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ABSTRACT 

The atmospheres of substellar objects contain clouds of oxides, iron, silicates, and other re- 
fractory condensates. Water clouds are expected in the coolest objects. The opacity of these 
'dust' clouds strongly affects both the atmospheric temperature-pressure profile and the emer- 
gent flux. Thus any attempt to model the spectra of these atmospheres must incorporate a 
cloud model. However the diversity of cloud models in atmospheric simulations is large and 
it is not always clear how the underlying physics of the various models compare. Likewise the 
observational consequences of different modelling approaches can be masked by other model 
differences, making objective comparisons challenging. In order to clarify the current state 
of the modelling approaches, this paper compares five different cloud models in two sets of 
tests. Test case 1 tests the dust cloud models for a prescribed L, L-T, and T-dwarf atmospheric 
(temperature T, pressure p, convective velocity V(;onv)-structures. Test case 2 compares com- 
plete model atmosphere results for given (effective temperature Ten, surface gravity log g). All 
models agree on the global cloud structure but differ in opacity -relevant details like grain size, 
amount of dust, dust and gas-phase composition. These models can loosely be grouped into 
high- and low-altitude cloud models whereas the first appear generally redder in near-infrared 
colours then the later. Comparisons of synthetic photometric fluxes translate into an modelling 
uncertainty in apparent magnitudes for our L-dwarf (T-dwarf) test case of 0.25 < Am < 0.875 
(0.1 < Am < 1.375) taking into account the 2MASS, the UKIRT WFCAM, the Spitzer IRAC, 
and VLT VISIR filters with UKIRT WFCAM being the most challenging for the models. 
Future developments will need closer links with laboratory astrophysics, and a consistent 
treatment of the cloud chemistry and turbulence. 

Key words: Stars: atmospheres - Stars: low-mass, brown dwarfs 



1 INTRODUCTION 

The atmospheres of L dwarfs are characterised by clouds, formed 
principally of silicate, oxide and iron grains, which shape their 
emergent spectra. Likewise the atmospheres of the early T dwarfs 
are distinguished by the progressive departure of cloud opacity. 
At even lower effective temperature, giant-gas planets are again 
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covered in clouds and other chemical components become impor- 
tant. Any attempt to derive fundamental properties of these objects 
from their spectra hinges on an understanding of the chemistry and 
physics of clouds. Yet clouds are inherently difficult to model since 
they can feedback into the chemistry and the physics of the entire 
atmosphere. Because of the complexity of this problem a number 
of independent groups have taken very different approaches to de- 
scribe the cloud formation and the cloud properties of substellar 
objects as a function of gravity, effective temperature, and metal- 
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licity. Here we make a first attempt to compare the quantitative pre- 
dictions of these various approaches in order to better understand 
the models themselves as well as the uncertainty which remains in 
application of these models to real objects. 

Atmospheric physics classically involves hydrodynamics, ra- 
diative and convective energy transport, and gas phase chemistry. 
Effects of magnetic fields are neglected. Ideally, the only free pa- 
rameters are the effective temperature Teff, the surface gravity g, 
radius R, or mass M,, and element abundances e;. In order to solve 
such a coupled system of equations in a computationally reasonable 
time, assumptions like the hydrostatic equilibrium, mixing length 
theory and chemical equilibrium are made. Inside substellar atmo- 
spheres, chemical equilibrium of the gas phase is justified due to 
high collision rates between gas-phase constituents. Irradiation or 
atmospheric flows may invalidate this assumption in the upper at- 
mospheric layers. The validity of hydrostatic equilibrium and mix- 
ing length theory have been studied in comparison to large eddy 
simulations for M-type stars (Ludwig et al. 2002, 2006) and we 
know from the direct observation of solar system giant planets at 
low gravities and effective temperatures that hydrostatic equilib- 
rium is an appropriate down to very low pressures (~ 1 yubar) in 
the atmosphere. The assumption of hydrostatic equilibrium cou- 
pled with the mixing length theory is computationally extremely 
efficient and an accurate approximation in particular if one is aim- 
ing at synthetic spectrum calculations. 

The striking difference of substellar atmosphere models com- 
pared to the classical stellar approach is the necessity to model 
the formation of clouds and their feedback onto the entire atmo- 
sphere. New physics needed to be considered and different tribes 
emerged being inspired by AGB star dust formation (Helling et al. 
2001a; Woitke & Helling 2003, 2004), by terrestrial cloud forma- 
tion (Ackerman & Marley 2001, Cooper et al. 2003), by measure- 
ments for solar system planets (Rossow 1978, Marley et al. 1999), 
or driven by practical considerations (Tsuji et al. 1996 a,b; Allard et 
al. 2001). The first attempts on cloud modelling in brown dwarf at- 
mospheres were undertaken by Lunine, Hubbard & Marley (1986) 
and Tsuji et al. (1996 a,b) who suggested the influence of clouds 
on the spectral appearance of brown dwarfs. See also Ackerman 
& Marley (2001) for a review and comparison of the earlier cloud 
literature. 

The overall, phenomenological understanding of cloud forma- 
tion in substellar objects has converged to the picture that dust (or 
condensates, see Table [3]( forms at a certain height in the atmo- 
sphere where it acts as an efficient element sink leaving behind a 
depleted gas phase. The departure of TiO and FeH spectral lines 
from the M to the L dwarfs testifies to this process. The dust then 
settles gravitationally taking condensed elements with it. Convec- 
tion and atmospheric mixing replenishes the condensing gas, re- 
sulting in a steady state. The details of this picture, however, leave 
room for debate. It is for example not clear where the dust pre- 
cisely starts to form since this depends on the details of the model 
assumptions. 

Current models generally employ one of two physical ap- 
proaches to understand this process. In the first paradigm gas is 
mixed upward into higher altitudes. Dust then forms, falls down 
and meanwhile grows until it evaporates below the cloud base 
(Woitke & Helling 2003). The second paradigm imagines the op- 
posite limiting case in which the gas is well mixed from the deep 
atmosphere only up to a cloud base. Grains and gas are transported 
above cloud base by mixing and grains fall down under the influ- 
ence of gravity (Ackerman & Marley 2001, Allard et al. 2003). 
These two branches rely on fundamentally different model assump- 



tions: (i) the phase-non-equilibrium concept of kinetic dust for- 
mation (Woitke & Helling 2003, 2004; Helling & Woitke 2006, 
Helling, Woitke & Thi 2008), and (ii) the phase-equilibrium con- 
cept of thermal stability (Tsuji et al 1996b, Tsuji 2005, Allard et 
al. 2001, Ackerman & Marley 2001, Burrows et al. 2006; Cooper 
et al. 2003). While (ii) represents the end-state which a dust form- 
ing system achieves for / — > oo, (i) describes the kinetic process 
of the formation of cloud particles on finite timescales limited by 
mixing and rain-out. The models also differ in the choice of dust 
materials which are assumed (i) to form or (ii) to be present in 
the atmosphere. Both areas need serious attention and correspond- 
ing material properties should be obtained either experimentally 
(see discussion Sect. |4j or from ab initio calculation (e.g. Jeong 
et al. 2000, Patzer et al. 2005) which both are beyond the scope 
of this paper. Given these model conceptions, a number of differ- 
ent model approaches have been developed to reproduce observed 
spectra (Tsuji et al 1996a, Tsuji 2005, Allard et al. 2001, Acker- 
man & Marley 2001, Burrows et al. 2006, Cooper et al. 2003, Dehn 
2007, Helling et al. 2008a,b) or providing detailed information on 
the dust complex itself (Woitke & Helling 2003, 2004; HeUing & 
Woitke 2006, Helling, Woitke & Thi 2008). 

Driven by this diversity in the field, the aim of this paper is 
to provide information and to perform comparative studies of mod- 
els that aim to describe the dust clouds in substellar atmospheres. 
Kleb & Wood (2004) demonstrated that such component-based test 
studies are an essential part of scientific methods. As the number, 
n, of model components increases, the interactions amongst them 
grow as 12. We therefore need to perform verifications on the 
components (here: cloud models) but also to use the method of 
manufactures solution^ (e.g. Kleb & Wood 2004) to verify that 
the entire system (here: model atmosphere) attains its theoretical 
order-of-accuracy properties. This goes beyond what has been and 
could be provided in the literature so far. 

Our paper begins with a summary of the various dust cloud 
models. We provide for the first time a comparative presentation of 
the different approaches concerning chemistry and dust modelling 
(Sect. |2). Based on a workshop held in Leiden in October 2006 

we present test cases where we first separate the components 
for chemistry and dust cloud modelling from the complete atmo- 
sphere problem (Sect. 13. It . This allows us to judge the order-of- 
accuracy properties of the complete models with respect to chem- 
istry and dust formation which both are essential ingredients for the 
solution of the radiative transfer problem. Section [331 demonstrates 
the results for the complete substellar atmosphere problem, syn- 
thetic photometric fluxes, and colours are calculated and synthetic 
trust ranges derived from independent models are given. 

Comparative studies have been carried out for simulations of 
radiative transfer (Pascucci et al. 2004, Iliev et al. 2006), of white 
dwarfs (Barstow et al. 2001) or photon dominated region (Rollig et 
al. 2007). No comparison study has been presented so far for sub- 
stellar atmospheres (brown dwarfs and planets). While we com- 
pare a number of model predictions, including emergent spectra, 
we refrain from comparing against spectra of individual substellar 
objects, since there are as yet no such objects with independently 



Manufactures solutions are tests with a known results. In our case, we 
manufacture the input quantities to a certain extent and compare the solu- 
tions (test case 1, Sect |3.1) . 

^ http://www.lorentzcenter.nl/lc/web/2006/203/info. php3?wsid=203 
' http://phoenix.hs.uni-hamburg.de/BrownDwarfsToPlanetsl/ 
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Table 1. Definition and units of quantities. The quantities plotted in Figs[T]- 
|9]are liigliliglited in boldface. 



Quantity 


Definition 


Unit 


T 


local gas temperature (Fig-Elll) 


K 




local gas pressure 


dyn cm"^ 


^conv 


convective velocity 


cm s"' 


e\ 


gas element abundance 


(i = H, He, ...) 





deep element abundance 




iS 


supersaturation ratio 




Pd/Pgas 


dust to gas mass ratio (Fig.[3 




Pi 


dust mass density 


gem"' 




= 2s ".s^s for a = const 










s 


dust species 






(e.g. Fe[s], Mg2Si04[s], . . .) 






number of dust particles of kind ^ 


-3 

cm 


a 


grain size 


cm 


Ms 


mass of dust particle of kind s 


g 


f(a) 


grain size distribution: 


cm^^'cm"' 




number of dust grains n per grain size 






a and per gas volume 




<a} 


mean particle size (Fig-ID 


cm 




Us X 












for fixed a = oq: 






f{a) = S(a - ao) =^ (a) = ao 




V /V, , 

s / tot 


volume fraction of dust kind s (Fis.[3i 


1 




total dust volume of dust kind s 


cm' 


tot 


total dust volume 


cm' 


F{A) 


surface flux (Fig.|6]|7) 


erg s"'cm"-A"' 


A 


wavelength 


/jm 


logFc 


photometric flux (Fig. [8) 


ergs cm 'A 




flux convolved with filter c 






f;'- F(A)tTanSc{A)dA 

JAi 






tans^(A)dA 




log Fco 


reference photometric flux (Vega) 




tranSc (A) 


function of a photometric filter c 






between Ai and A2 (see Table[6) 




mi — m2 


colour (Fig.|9) 






fc2(Ai2) , fcl(Ail)^ 

= 2.5 log log 






^ ^F,02(Ai2) ^F.oi(A^i)^ 




m 


apparent magnitude 






= -2.5 log X'|- F(i)transc dA 




constrained mass, age, and metallicity against which spectral mod- 


els can be compared. 





2 APPROACHES TO CHEMISTRY AND DUST CLOUD 
MODELLING IN BROWN DWARF ATMOSPHERES 

2.1 Gas phase chemistry 

Each of the cloud models to be summarised (Sect. 12.2b assumes 
local thermodynamic equilibrium (LTE) when modelling the gas 
phase chemistry. In stellar atmospheres, departures from LTE can 
arise from interactions of atoms and molecules with the non- 
thermal radiation field (Woitke, Kriiger & Sedlmayr 1996) but this 
effect is negligible in dense substellar atmospheres (Hauschildt et 
al. 1997; Schweitzer, Hauschildt & Baron 2000). Dust and gas are 



assumed to have the same temperature {Ti^st = T^gas)- The gas phase 
abundances determine the kind and the amount of dust condensing 
and are in turn determined by the amount of elements not bound by 
the dust. All codes use equilibrium constants Kp in their gas-phase 
treatments (Fegley & Lodders 1994; Tsuji et al. 1996b, Tsuji 2005; 
AUard et al. 2001, 2003; Lodders 2003; Helling & Woitke 2006; 
Helling, Woitke & Thi 2008). In reality, differences may arise due 
to the selection of the input data, which would be apparent in direct 
comparisons among the various gas-phase models. 

The aim of this paper is to investigate and quantify the dif- 
ferences arising from different cloud model approaches rather than 
testing thermodynamic data sources. While we do not expect large 
uncertainties due to possible difference in thermodynamical gas- 
phase data, the results of the gas-phase chemistries used by differ- 
ent modellers will differ if different sources for element abundances 
were used (see discussion in Sect. |3.U . 



2.2 Dust cloud models 

In the following we summarise five different cloud models which 
are used in substellar atmosphere simulations, and which are in- 
volved in our comparative calculations (Sect. [T2ll3.4t . While there 
are many differences, ultimately all of the models face the same un- 
derlying physical challenges. We will try to note conceptual simi- 
larities and differences as we describe the models below. 
The descriptions include where appropriate: 

• the link between the cloud module and the atmosphere code, 

• the physical ideas and their representation, 

• the treatment of the cloud chemistry. 



2.2.1 Tsuji model 

2.2.1.1 Dust and gas-phase treatment in model atmosphere 
code: Tsuji and collaborators apply the methods of non-grey ra- 
diative transfer (in hydrostatic and radiative-convective equilibrium 
under LTE) to dusty photospheres with almost no modification, ex- 
cept that the solid and liquid phases are considered in addition to 
the gas phase in chemical equilibrium. In solving chemical equi- 
librium, the Tsuji models aim to provide the abundances of ions, 
atoms, molecules, and dust grains that contribute to the opacities 
rather than to derive a complete solution for all elements. Thirty- 
four elements are considered in charge conservation, 16 elements 
(H, C, N, O, Na, Mg, Al, Si, P, S, CI, K, Ca, Ti, V, Fe) in molec- 
ular formation, and 8 elements (Mg, Al, Si, Fe, K, Ca, Ti, V) in 
dust formation. The chemical equilibrium computation includes 
83 molecules and is based on a previous examination of about 
500 molecular species (as for details, including the thermochem- 
ical data, see Tsuji 1973). Dust grains composed of Fe, Si, Mg, 
and Al in form of metallic iron, enstatite (MgSiOs), and corundum 
(AI2O3) are considered as sources of dust opacity. The abundances 
are solved as being in phase equilibrium with gaseous species. No 
other dust species composed of Fe, Si, Mg, or Al are considered 
for simplicity. Also, abundances of some gaseous species impor- 
tant as sources of gaseous opacity suffer large reduction by the 
dust formation, and such effects are approximated by perovskite 
(CaTiOj) for Ti, melilite (CajMgSizOv) for Ca, VO(cr) for V, and 
K2S(cr)/KCl(cr) for K, since gaseous TiO, CaH, VO, and K are 
important sources of gas opacity. 
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Table 2. Dust cloud models in substellar atmospheres (z - atmospheric height; s - dust species). 

The references are the following: (S) - Tsuji (2000); ® - Tsuji et al. (1996b); ® - Tsuji (2002, 2005), ® - Allard et al. (2001); ® - AUard et al. (2003); ® - 
Ackerman & Marley (2001), ® - Woitke & Heffing (2003), Helling, Woitke & Thi (2008) 

Author Assumptions Model variants 



grain size a grain composition supersaturation references 



Tsuji a = lO'^jum homogeneous S = 1 case B full dusty model (D 

case C dust cleared model @ 

UCM dust between ® 

AUard & Homeier f(a) = a^^-^ (ISM) homogeneous S = 1 dusty full dusty model ® 

cond dust cleared model ® 

time scale comparison homogeneous 5 = 1.001 settl ® 

(/ Rossow 1978) 

Ackerman & Marley log-normal /(a, z) homogeneous S = 1 sedimentation ® 

efficiency 

Helling & Woitke f{a,z) dirty S=S{z,s) ® 



2.2.1.2 Cloud model: The Tsuji models assume that dust forms 
in the photosphere as soon as the thermodynamical condition for 
condensation is met, i.e. the supersaturation 5 = 1. Then the lay- 
ers cooler than the condensation temperature (Tcond) are assumed 
to be filled with dust grains (case B) which act as element sink 
and opacity source. Another extreme case is that the dust grains all 
precipitate as soon as they are formed and the atmosphere is thus 
clear of dust {case Q, hence the dust acts as element sink but not as 
opacity source. Finally an intermediate case (the "Unified" or UCM 
case) in which grains condense at T^oni, but precipitate at a slightly 
lower temperature termed the critical temperature, T„ is also con- 
sidered. In this case the dust cloud appears in a restricted region of 

Tcr < Z" < T^cond. 



While Tcond is well defined by thermal stability, is left as a 
free parameter to be estimated empirically. If T^r is equal to Tcond, 
all the dust grains will precipitate as soon as they are formed {case 
C). On the other hand, if is as low as the surface temperature, all 
the dust grains formed will survive in the fully dusty photosphere 
{case B). If Jcr differs only slightly from Tcond, the dust cloud will 
be quite thin while the dust cloud will be rather thick if T^r is much 
lower than Tcond. Thus Ta is essentially a measure of the thickness 
of the dust cloud and thus has a significant effect on the infrared 
colours (not unhke fsed in the Ackerman & Marley (2001) model). 
As a free parameter (along with Tefs) can be inferred from the 
observed infrared colours. For this purpose, reasonably accurate 
values of Tf.g can be inferred from the luminosities based on the 
observed parallaxes and bolometric fluxes (e.g. GoUmowski et al. 
2004; Vrba et al. 2004). But it appeared that the infrared colours 
differ significantly even for the same T^^^ (e.g. Marley et al. 2005; 
Tsuji 2005) and this fact implies that T^, also differs for the same 
Tf.g. For example, four cool dwarfs from spectral type L6.5 to T3.5, 
whose infrared spectra are quite different, appear to have almost the 
same empirical Jeff at about 1400+100 K. Such very different spec- 
tra of almost the same T^ff could be explained reasonably well by 
assuming different values of Tc„ i.e. dififerent thickness of the dust 
cloud (see Fig. 10 of Tsuji 2005). 



In the Tsuji models all grains have radius a = 0.01 /im. In 
the Umit of such small sizes the dust opacity is independent of the 
particle size for a fixed mass of dust. 



2.2.2 Allard &■ Homeier model 

2.2.2.1 Dust and gas-phase treatment in model atmosphere 
code: The Allard & Homeier models solve for chemical equi- 
librium in the gas phase by minimisation of the functional errors, 

where the functions are the elemental and charge conservation, 
Saha equation, and mass-action law for each of 40 elements, with 
up to 5 ionisation levels per atom, and for some 600 molecules and 
nearly as many condensate using thermochemical data from many 
sources including a compilation of the JANAF tables (Chase et al. 
1986; for details see Allard et al. 2001 and Allard & Hauschildt 
1995). 

Allard et al. (2001) modelled the limiting effects of cloud for- 
mation {dusty, cond) on the spectral properties of late M and L to T 
brown dwarfs by treating dust in chemical equilibrium with the gas 
phase. For the grains construction and opacities in the dusty and 
the cond models, an interstellar size distribution of spherical and 
chemical homogeneous grains is assumed. A slight supersaturation 
(1.001) is assumed in the setd models (Allard et al. 2003, 2007). 

The grain sizes are calculated from the comparisons between 
time-scales for mixing due to convective overshooting as prescribe 
by (Ludwig et al. 2002) and condensation and gravitational settling 
according to Rossow (1978). The thermal structure of the model at- 
mosphere is solved on a fixed optical depth grid at 1.2/im assuming 
no dust opacity contribution. The dusty models accounted for dust 
opacity while the cond models did not (Allard et al. 2001). The settl 
models involve a detailed cloud model (Allard et al. 2003, 2007) 
which is solved for the thermal structure of the atmosphere to find 
the grain size and abundances distributions as function of depth. 
In order to account for the gas cooling effects as it is propelled by 
convective turbulence from the top of the convection zone towards 
the top of the atmosphere, the cloud model is solved by depleting 
gas phase abundances layer-by-layer from the innermost (assumed 
of solar composition) to the outermost layer. 

The resulting stratified elemental abundances and number 
densities of species are then used in the radiative transfer solver ap- 
plying the Mie equation and complex refractive index for calculat- 
ing the dust opacities (Ferguson et al. 2005). Models are converged 
removing thereby any possible cloud opacity inconsistencies be- 
tween thermal structure and radiative transfer. 

2.2.2.2 Cloud model: For the settl model atmosphere (Allard et 
al. 2003, 2007) in each layer, the condensation, sedimentation and 
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Table 3. A brief dictionary for multiple meanings and phrases (a, b, . . .) used by different authors 

(T - Tsuji; AH - Allard & Homeier; MAL - Marley, Ackerman & Lodders; HW - Helling & Woitke; R - Rietmeijer). 



dust 


a) 


(HW) 


general term for small solid particles, grains, liquid droplets 




b) 


(MAL) 


condensate 


condensation 


a) 


(T, AH) 


dust formation by conversion of vapour to solid (or liquid) particles 


nucleation 


a) 


(MAL, HW) 


seed particle formation 


growth 


a) 


(HW) 


formation of condensate species by chemical surface reaction on an existing surface 


evaporation 


a) 


(^rlVV ) 


reverse growth process (Tgvap ■ evaporation time scale) 


drift 


a) 


(HW) 


relative motion between gas and dust 




b) 


(MAL) 


gravitational settling 




c) 


(MAL) 


rain, rain-out 




d) 


(AH, MAL) 


sedimentation: the falling of cloud particles under the influence of gravity (Tsej - sedimentation time scale) 




e) 


(T, MAL) 


precipitation: formation of cloud particle for which Tevap s> Tjed (from Rossow 1978) 


homogeneous 


a) 


(HW) 


seed formation by addition of the same monomer species forming clusters of increasing sizes until they 


nucleation 






achieve solid state character 




b) 




formation of first condensate that will grow to increasingly larger clusters 


heterogeneous 


a) 




seed formation by addition of different monomer species 


nucleation 


b) 


(R) 


formation of condensed species onto an existing seed (y growth) 


primary condensate 


a) 


(MAL) 


a condensate forming from the gas by gas-gas reactions 


secondary condensate 


a) 


(MAL) 


a condensate forming by gas-solid reaction with previously existing solid or liquid phase 


coagulation 


a) 




formation of one particle from two colliding cloud particles 


coalescence 


a) 


(AH) 


coagulation caused by size dependent sedimentation velocities of cloud particles of different sizes 


element conservation 






exchangeable used with mass balance since no elements should be created inside the atmosphere 


supersaturation 


a) 


(HW, T, AH) 


S , the ratio of the monomer particle pressure to the saturation vapour pressure 








(see Helling, Woitke & Thi 2006, Appendix for discussion) 




b) 


(MAL) 


S - 1 , the vapour pressure in excess of saturation divided by the saturation vapour pressure 



coalescence timescales (Rossow 1978; see Table[3]l are compared to 
the mixing timescale prescribed by Ludwig et al. (2002) as follows: 

a) the equilibrium size between mixing and sedimentation is cal- 
culated and the growth time scale (condensation and coalescence) 
is computed for that size; 

b) the mixing time scale is then compared to the growth time 
scale: 

• if growth is faster the condensates are found to be depleted, 
and the fraction of condensates is recomputed so as to obtain a 
growth time scale equal to the mixing time scale; 

• when mixing is faster, the growth is limited by the replen- 
ishment with fresh condensable material from deeper layers, and, 
while the condensate fraction is stable, a mean size is recom- 
puted corresponding to an equilibrium between mixing and con- 
densation. 

Given the new cloud description the elemental abundances 
are then readjusted which produces a new equilibrium condensate 
fraction. These steps are repeated until the condensate fraction no 
longer changes. This is a time consuming process which however 
guaranties that the chemistry reflects the cooling path of the gas. 

Another essential input to the model is the description of the 
mixing timescale. Within the lower classically convective unstable 
atmosphere layers a mixing velocity is readily obtained from the 
results of mixing length theory, which is implemented in PHOENIX 
in the formulation of Mihalas (1978). Since the principal cloud for- 
mation region is located well above the Schwarzschild boundary, 
one is confronted with the task of extrapolating the velocities over 
several pressure scale heights. For a phenomenological description 
of this velocity field this group draws on the results of the hydro- 
dynamical simulations of late M dwarfs by Ludwig et al. (2002), 
which show in general an exponential decline of mass transfer by 



overshooting with decreasing pressure, after an initial transition 
zone. Further simulations by Ludwig (2003) indicate a steepen- 
ing of this decline with surface gravity. The mass exchange fre- 
quency following these simulations is parameterised in analogy to 
the Helling & Woitke -model as log T,nix(z) = log Tmlt + y6(log Po - 
log p(z), where the base value of the mixing time scale (where 
p = Po) within the convectively unstable layer is given by its mix- 
ing length theory value Tmlt = oHp/Vcom and a = 2.0 the mixing 
length parameter (Ludwig et al. 2002). The slope can be derived 
from Ludwig (2003) as p = l^Jg^, where g5 is the surface gravity 
in units of 10^ cm s"^. Since the calibration of this relation involved 
an extrapolation of the M dwarf simulations to lower Teff , the mod- 
els allow for adjusting the slope /? by a factor of up to 3, adopting a 
factor of 1 for the 1800K and 1400K test cases, and 2 for the 600K 
and lOOOK cases. 



2.2.3 Marley, Ackerman &■ Lodders model 

2.2.3.1 Dust treatment in model atmosphere code: The Mar- 
ley, Ackerman & Lodders - modelling treats the upward convective 
mixing of a gas, its condensation, and the sedimentation of con- 
densate through the atmosphere of an ultracool dwarf. The compo- 
sition and cloud structure at each point in a trial atmosphere model 
is computed, based on the existing profile and then this informa- 
tion is used to iterate towards the next trial profile. The chemistry 
at each pressure/temperature point is interpolated within a table of 
atomic and molecular abundances computed for chemical equilib- 
rium (Freedman, Marley & Lodders 2008). The cloud is computed 
by applying the Ackerman & Marley (2001) cloud model. 

2.2.3.2 Gas-phase treatment: Abundances of gas species are 
calculated with the CONDOR code (Lodders & Fegley 1993; Fe- 
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gley & Lodders 1994; Lodders 2003) which calculates chemical 
equilibrium compositions by considering the dual constraints of 
mass balance and chemical equilibrium. Input data required for 
the code are thermodynamic properties of the gas-phase species 
and compounds (e.g., equilibrium constants), appropriate elemen- 
tal abundance tables for the system, temperature and total pressure. 
The equilibrium constants used in the CONDOR code are com- 
puted from the Gibbs free energy (AG) data, which are directly 
proportional to the logarithm of the equilibrium constant (In K^) as 
AG = -RT In A'p {T - gas temperature, R ~ gas constant). The code 
considers ~ 2000 gas species (including ions) and ~ 1700 solids 
and liquids for compounds of all naturally occurring element^ 

2.2.3.3 Condensate treatment: For application to substellar at- 
mospheres, the CONDOR code treats condensate formation by re- 
moving primary condensates (i.e., condensates that form from con- 
densing gases) from the gas into cloud layers (Lodders 2004, Lod- 
ders & Fegley 2006; see Table[3}. An important consequence of this 
approach is that secondary condensates arising from gas-solid re- 
actions as would be predicted by pure equilibrium are excluded be- 
cause the primary condensates are assumed to settle into clouds and 
are thus no longer available for reaction with the cooler gas above 
the clouds. For example the computation assumes that iron grains (a 
primary condensate) do not react with H2S gas to form FeS at lower 
temperatures (~ 700 K) where the secondary FeS would form if Fe 
metal were still present. Instead the H2S remains in the gas phase as 
is observed in the deep atmosphere of Jupiter (Niemann et al. 1998) 
where H2S is only removed into NH4SH clouds below ~ 200 K (Fe- 
gley & Lodders 1994; Visscher et al. 2006). Likewise a detection 
of H2S in a cool T dwarf would confirm the inhibition of secondary 
condensation. Marley et al. (2002) argue that the far red colours 
of T dwarfs can only be reproduced if secondary condensation of 
alkali-bearing phases is indeed inhibited. 

Both the cloud model of the settling of primary condensates 
and the chemical equilibrium model assume that at a given tem- 
perature below the condensation temperature the gas phase abun- 
dances of the elements sequestered by condensation are set by the 
respective vapour pressure of the primary condensate. With this 
common assumption the cloud and chemical computations are fully 
self-consistent. 

2.2.3.4 Cloud model: The cloud model (Ackerman & Marley 
2001) parameterises the efficiency of sedimentation of cloud parti- 
cles relative to turbulent mixing through a scaling factor, f^^ (Eq. 
4 in Ackerman & Marley 2001). Large values of /sed correspond 
to rapid particle growth and large mean particle sizes. In this case 
sedimentation is efficient, which leads to physically and optically 
thin clouds. When f^^^ is small particles are assumed to grow more 
slowly and the amount of condensed matter in the atmospheric is 
larger and clouds thicker. In this sense small /sed is somewhat com- 
parable to the Tsuji models with a large difference between Tcond 
and Tcii while large /sed is similar to the opposite case. Unlike the 
Tsuji models Ackerman & Marley (2001) compute a particle size 
profile for each condensate in each model atmosphere. 

For a fixed atmospheric profile, /ed and a description of the 
width of the particle size distribution, the Ackerman Sl Marley 
(2001) model uniquely predicts the variation in mean particle size 
and particle number density through the atmosphere. Thus families 

All elements up to Bi (number 83) excepting Tc and Pm but adding Th 
and U, for a total of 83. 



of models, i.e. set of models with varying sets of parameters Teff, 
logg, fsed, each with a unique /ed, can be produced. No attempt 
is made to model microphysical processes of dust growth and co- 
agulation. Instead it is assumed that the micro-physical processes 
acting within the cloud are able to produce the particle sizes im- 
plied by any specified value of /ed- 

In terrestrial rain clouds the particle size distribution is often 
double-peaked (Ackerman et al. 2001), with small particles that 
grow from condensation of the vapour co-existing with larger drops 
that have grown by by collisions between particles. In the Acker- 
man& Marley (2001) approach a single, broad log-normal particle 
size distribution is intended to capture the likely existence of such 
a double-peaked size distribution. They employ a width cr = 2 for 
all cases, although this can in principle be varied. 

Like the Allard & Homeier - model, the Marley, Ackerman 
& Lodders cloud - model use the mixing length theory to compute 
gas velocities in the convection zone and must employ some other 
description to specify velocities above the radiative-convective 
boundary. They describe mixing in radiatively stable layers by 
specifying an eddy diffusion coefficient, A'eddy = H'/Tmix, where H 
and T„,ix are the scale height and mixing time. Experience with the 
radiative stratospheres in the solar system (e.g., Atreya et al. 1991; 
Bishop et al. 1995; Moses et al. 2004) shows that typical values of 
A'eddy in these atmospheres lie in the range of lO** to 10^ cm^s. Ob- 
servations of ammonia and CO in the atmosphere of the T7.5 dwarf 
Gliese 570D imply K^My ~ cm^s (Saumon et al. 2006; Geballe 
et al. 2008). Comparisons of the mid-infrared colours of L dwarfs to 
models that include chemical mixing (Leggett et al. 2006), suggest 
^cddy ~ 10* cm' s"'. The Marley, Ackerman & Lodders - models 
reported here set A'eddy > 10^ cm- s"' at all points in the atmosphere 
with a smooth transition from the convective zone to this value. 

As with variation of cloud thickness to match variation in J-K 
at fixed Teff in the Tsuji models, changes in /ed produce atmosphere 
models with a range of near-infrared spectra and colours. Burgasser 
et al. (2007) and Stephens et al. (2008) have shown that the spec- 
tra of bluer-than-average L dwarfs can be fit by models employing 
large /ed while redder-than-average L dwarfs seem to require small 
/ed. The Marley, Ackerman & Lodders reproduce spectra across 
the L to T transition by employing models with progressively larger 
/ed with later spectral type (Gushing et al. 2008; Stephens et al. 
2008; Leggett et al. 2008). 

2.2.4 Helling &■ Woitke model 

I.IAA Dust treatment for a model atmosphere code: The 

Helling & Woitke approach is fundamentally different from the 
previous models in two important ways. First, this model follows 
the trajectory of an ensemble of dust grains downwards from the 
top of the atmosphere instead of upwards from the bottom. This 
approach is based on the phenomenological analogy to thunder- 
storm where large air masses are advected upwards before rain- 
drops do form. Dust clouds in substellar objects are considered 
stationary, i.e. uncondensed gas is mixed upward from which dust 
particles continuously form, settle gravitationally, and evaporate. 
In this stationary situation, the downward directed element transfer 
via precipitating dust grains is balanced by an upward mixing from 
the deep interior by convective and overshoot-motions (Helling et 
al. 2001b, Helling 2003). The second major difference from the 
other approach is that the Helling & Woitke approach kinetically 
describes the cloud particle formation as phase-transition process 
by modelling seed formation, grain growth/evaporation, sedimen- 
tation in phase-non-equilibrium, element depletion, and their in- 
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teractions. Dust moment equations describing tliese processes are 
derived from rate equations and are solved as a function of fieight 
z for a given (r, p, iVonv) atmospliere structure. Tlie equations are 
integrated inward. 

I.IA.I Gas-phase treatment: The composition of the gas 
phase is calculated assuming chemical equilibrium for 14 elements 
(H, He, C, N, O, Si, Mg, Al, Fe, S, Na, K, Ti, Ca) and 158 molecules 
with equilibrium constants fitted to the thermodynamical molecu- 
lar data of the electronic version of the JANAF tables (Chase et al. 
1986). The equilibrium constant for TiC are from Gauger et al. (see 
Helling et al. 2000), for CaH from Tsuji (1973), and FeH from Bur- 
rows (priv. com.). Solar elemental abundances are assumed at the 
lower boundary of the model atmosphere, assuming a well-mixed 
gas-phase solar composition, and first ionisation states of the ele- 
ments are calculated. Element conservation equations are auxiliary 
conditions which take into account the loss of elements in the gas 
phase by nucleation, growth, and drift and the gain by evaporation 
(Woitke & Helling 2004). 

2.2.4.3 Condensate treatment: The condensates considered 
during the solution of the Helling & Woitke dust model equations 
are treated in full phase-non equilibrium. The supersaturation ratio 
S is calculated from the gas phase composition in chemical equi- 
librium. 

I.IAA Cloud model: The dust formation starts with the for- 
mation of seed particles (nucleation). The nucleation rate is cal- 
culated for homogeneous (Ti02)A/-clusters applying the modified 
classical nucleation theory (Gail et al. 1984; see Eq. 34 in Helling & 
Woitke 2006). The calculation of the nucleation rate relies on quan- 
tum mechanical calculations for the formation of Ti02-seeds by a 
step-wise addition of Ti02 molecules (Jeong et al. 2000). The nu- 
cleation rate determines the number of dust particles. These seeds 
grow to macroscopic sizes by gas-solid surface reactions. Because 
many compounds can be thermally stable almost simultaneously 
in substellar atmospheres, the simultaneous growth of 12 solids 
Ti02[s], SiO[s], SiOjM, Fe[s], FeO[s], Fe203[s], FeS[s], MgO[s], 
MgSi03[s], Mg2Si04[s], Al203[s] and CaTiOjIs] by 60 surface re- 
actions onto Ti02-seed particles is modelled (Helling, Woitke & 
Thi 2008). These dirty grains are modelled to be composed of a 
homogeneous mix of numerous islands of the different, pure con- 
densates (Helling & Woitke 2006). Drift transports existing par- 
ticles into region where they might continue to grow before they 
evaporate in the deeper, warmer atmosphere. While reactions on an 
existing grain surface proceed if the gas is supersaturated (5 > 1) 
with respect to this particular reaction (Helling & Woitke 2006), 
the seed formation can only take place when the gas is highly su- 
persaturated (5 » 1). If the gas is under-saturated (5 < 1) the solid 
will evaporate. 

The majority of dust grains that build up the cloud layer are 
found in a subsonic gas for which Knudsen numbers are small 
(see Woitke & Helling 2003). The respective kinetic description is 
solved in form of conservation equations which allows a simultane- 
ous treatment of nucleation, growth, evaporation, drift, and element 
replenishment. The dust formation is modelled by applying conser- 
vation equations of dust moments Lj = J VJ''/(V)dV with /(V) 
the grain size distribution function. Nucleation, growth/evaporation 
and gravitational settling are source terms of these equations 
(Woitke & Helling 2003, Helling & Woitke 2006). The solution 
of the dust moment equations and element conservations deter- 
mines quantities like grain sizes, grain material composition, total 



grain volume, remaining gas-phase element abundances. The ele- 
ment replenishment is treated by introducing a parameterised mix- 
ing time scale Tn,jx(z). Ludwig et al. (2002, 2006) show that, gen- 
erally speaking, the convectively excited hydrodynamical motions 
- and thereby the mixing - decay exponentially with increasing 
height above the convectively unstable zone resulting in an expo- 
nential decrease of the mass exchange frequency in the radiative 
zone from which Tmix(z) is derived as log Tmix(z) = log t™"(z) + P ■ 
[ 0, log Pa - log p{z) I with pq the pressure at the upper edge of the 
convective unstable zone, t"™"(z) = a/i/pVconv (q- = 2.0) the mini- 
mum value of the mixing time-scale occurring in the convectively 
unstable region and yS = A log /exchange/^ log p a 2.2. 

2.2.5 Dehn &■ Hauschildt + Helling &■ Woitke model 

2.2.5.1 Dust treatment in model atmosphere: The dust cloud 
model of Helling & Woitke (Sect. 12.2.4b has been adopted as 
module in the static Phoenix model atmosphere code (Dehn 
2007; Helling et al. 2008 a,b). the dust module receives the 
(r(z),p(z), Vconv(z)) structure from Phoenix and provides the dust 
number density, the solid's volume fractions, the mean grain size, 
and the remaining element abundances in the gas phase for each at- 
mospheric layer. Effective medium and Mie theory are then used to 
calculate the dust opacity in addition to the usual gas-phase opac- 
ity calculations. The temperature structure is found iteratively by 
a modified Unsold-Lucy correction algorithm. The adjusted atmo- 
sphere structure, including the solution of mixing-length theory to 
find Vconv> is an input for the dust module in the next iteration. Com- 
pared to the classical Phoenix solution, the computing time has in- 
creases by a considerable factor of since for each temperature iter- 
ation the dust module is called. The dust module itself iterates to 
solve the dust moment equation by fulfilling the element conserva- 
tion auxiliary condition. 

2.2.5.2 Gas-phase treatment, condensate treatment, cloud 
model: The gas-phase composition is calculated assuming chemi- 
cal equilibrium as described in Sect. |2.2.2l The cloud model is a re- 
duced version of Sect. 12.2.41 in order to keep the computation time 
effortable: The simultaneous growth of 7 solids Ti02[s], Si02[s], 
Fe[s], MgO[s], Mg2Si03[s], Mg2Si04[s], Al203[s] onto the TiOj- 
seed particles is considered. Only 32 surface reaction are taken into 
account (Dehn 2007). 



3 TEST CASES 

We explored the characteristics and capabilities of the cloud models 
summarised in Sect. [2] with two sets of test calculations. The test 
case 1 is the component-based test study (Kleb & Wood 2004), 
and it is designed to compare the dust cloud models alone by sep- 
arating them from hydrodynamics and radiative transfer treatments 
(including opacity calculations). Test case 2 utilises the method of 
manufactured solutions (Kleb & Wood 2004), and compares the 
results of completely iterated substellar model atmosphere simula- 
tions. Table[T]contains the definitions of the quantities discussed in 
the following. 

3.1 Test case 1: local quantities given 

Each dust cloud model is calculated for a prescribed set of 
(T, Pgns, , Vconv) or (r, pg^, ef, Fcom) profiles with T the local gas 
temperature, pgas the local gas pressure, v^om and Fconv being the 
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Figure 1. Input model structure (T, /.), i^onv) for Tefj =1800, 1400, 600 K 
with log g=5.0 and solar element abundances. 



convective velocity and the convective flux, respectively (Fig. [Hi. 
The deep, well-mixed element abundances e!* (i=H, Si, Mg, Ti, . . .) 
in the inner atmosphere have been chosen as solar according to 
Grevesse, Noels & Sauval (1992). 

Remarks on the solar element abundances: The solar element 
abundances published by Anders & Grevesse (1989), Grevesse, 
Noels & Sauval (1992), and Grevesse & Sauval (1998) are under- 
going considerable revisions but an agreement on the most correct 
value according to present knowledge has not yet been reached. 
The oxygen element abundance, as the most prominent example, 
has been revised downward to 6^ = 8.66 ± 0.05 by Asplund et 
al (2004) and to only = 8.76 ± 0.07 by Caffau et al (2008), both 
based on 3D hydrodynamical simulations of the solar photosphere 
in combination with non-LTE line transfer (see discussion in Caffau 
et al. 2008). However, Ayres et al. (2006) suggest ej^ = 8.84 from 
their measured solar CO lines. Additionally, the downward revision 
of the oxygen abundances greatly increases the difference between 
the internal sound speed predicted by solar models and the sound 
speed inferred from helioseismology (Christensen-Dalsgaard et al. 
2008). The determination of solar element abundances is a funda- 
mental problem for atmosphere physics and chemistry, and the final 
amount of dust formed in a cloud will depend on the element abun- 
dance values. However, the test of its implications goes beyond the 
scope of this paper. Therefore, test case 1 (Sect. |3.2| l applies the 
Grevesse, Noels & Sauval (1992). The element abundances used 
for test case 2 (Sect. l34t are listed in Table[5]for each of the atmo- 
sphere codes. 



Remarks on the mdcing time-scale t^\^: A mixing time-scale t^\,^ 
enters all cloud models except the Tsuji-model. Each model (Mar- 
ley, Ackerman & Lodders, Allard & Homeier, Helling & Woitke) 
does interpret, and hence, uses T„,ix in different ways. In the Mar- 
ley, Ackerman & Lodders-model, the sedimentation efficiency is 
parameterised relative to a mixing time-scale through a scaling fac- 



tor /sed, the Allard & Homeier-model includes a t^ix in a time-scale 
comparison to determine local mean grain sizes, and a r^ix influ- 
ences the rate of seed formation, the growth and the settling pro- 
cess since it enters a set of conservation equations in the Helling & 
Woitke-model. Hence, we refrain from directly comparing the mix- 
ing time-scales of the different cloud models, since a comparison of 
Tmix would be rather misleading regarding the cloud properties de- 
rived by each of the modeller groups. 

We have chosen to compare the models for the following stellar 
parameter which can be considered as examples for the L-, L-T, 
and T~dwarf atmospheres: 

L - dwarf Teff = 1 800K (provided by M. Dehn) 
L-T dwarf Teff = 1400K (provided by M. Dehn) 
T - dwarf Teff = 600K (provided by M. Marley) 

All models have logg = 5.0. The given (T,p) and (Vco„v,p) 
structures are shown in Fig. [T] The Teff = 600K model is consider- 
ably less extended in logp than the hotter models. Its convective 
velocity is very small. Therefore, a much less efficient convective 
overshooting is anticipated, and (for those models which assume 
gaseous transport) a less efficient element replenishment of the up- 
per atmospheric layers. 

3.2 Results: test case 1 

We compare four essential results of our dust cloud models which 
are needed for the opacity calculations in a complete atmosphere 
model: 

• dust content 

• mean particle size 

• dust material composition 

• gas-phase composition 

3.2.1 Dust content in the atmosphere 

We measure the dust content in the atmosphere by the dust-to-gas 
mass ratio Pdust/Pgas (for definition see Table[TJ Fig. |2]left). How- 
ever, the phase-equilibrium models allow two interpretation of this 
quantity, namely, the amount of dust acting as opacity source and 
the amount of dust acting as element sink (compare column 6 in 
Table Is}- Figure [2] depicts Pdust/Psas for dust opacity sources, and 
demonstrate for the Allard & Homeier-model that the difference to 
Pdus(/Pgas for dust element sinks can be significant. 

The dust-to-gas mass ratio shows where most of the dust is 
located in the cloud layers, and the extension of the cloud layer(s) 
differs for all models. The innermost cloud layers generally con- 
tain the maximum amount of dust, except in the Tsuji results where 
the clouds extension varies with T^r. All models have the same lo- 
cation of the inner boundary of the cloud (cloud base), since it is 
determined mainly by thermal stability. 

The maximum Pdust/Pgas is of about the same order of magni- 
tude for all models but the exact values differ (see Table |4j. Note 
that the amount of dust entering the radiative transfer calculation is 
usually smaller than the amount of dust causing the gas phase de- 
pletion in phase-equilibrium models (compare Allard & Homeier- 
model: gray lines in Fig. [2] left). The maximum Pdust/Pgas value 
is reached at different atmospheric altitudes in the different mod- 
els and it retains its value over different atmospheric extension. 
The Tsuji-models suggest the highest amount of dust in the atmo- 
spheres for Tcy = 1700K. The Allard & Homeier-models suggest 
the lowest amount of dust. The Teff = 600K test case is challenging 
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Figure 2. Test case 1 results for prescribed atmospheric structures for Teff =1800, 1400, 600K, all log g = 5.0 and solar element abundance. Left: Dust-gas-ratio 
pd/Pn Right: Mean particle size (a) [jim] 

Note: For Tsuji two cases of Ta are plotted for each Teff ; = 1900K - light colours (orange/light blue), T^^ = 1700K - dark colours (red/blue). For Marley, 
Ackerman & Lodders Pd/Pg = (PA1203 + PFt + PMg2Si04 + PH2o)/Pg (/sed = 2). The homogeneous H2O-, AI2O3-, Mg2Si04-, and Fe-particle have different 
sizes (different line styles). For AUard & Homeier Pd/Pe includes all species contributing to the depletion of the gas phase. The gray lines show their values 
which enters the radiative transfer calculation. For the Allard & Homeier Tejf = 600K-model, the opacity-species-only-log Pd/Pg values fall below the axis 
range depicted. For Helling & Woitke, the code has difficulties calculating clouds for T^fj = 600K in the inner atmosphere. 



for all models: No data could be provided from the Tsuji-model, 
the Allard & Homeier Pdust/Pgas have a local minimum, and the 
Helling & Woitke-model reaches the shallowest depth. The Mar- 
ley, Ackerman & Lodders Teff = 600K-model (solid green line in 
Fig.|2]l shows two well-separated cloud layers: a water layer high 
up in the atmosphere at a 10 dyn/cm^ and a silicate layer between 
10^ ... 10^ dyn/cm^ (compare Sect. I3.2."3] l. All other models pro- 
duce only the silicate layer. 



3.2.2 Mean particles sizes in the cloud layer 

Figure|2](right) shows the results for the means grain sizes {a) (def- 
inition see Table [TJ calculated for given {T, p, v^-am) profiles. The 
mean grain sizes are different amongst all models which reflects 
the different model assumptions made. Also, the grain size distri- 



bution function f(a) used to determine (a) is different in each of 
the dust cloud models (see Table[2ll. 

A common feature for all models is that small mean particle 
sizes (a) < 10^^/jm populate the upper cloud regions, except in 
the Allard & Homeier-model. This small grain size in the upper 
cloud layers are associated with very small dust-to-gas ratios of 
Pdust/Pgas < 10"*. Particle sizes increase inward and reach a certain 
maximum size (Marley, Ackerman & Lodders; Helling & Woitke), 
or are constant by assumption in the entire cloud (Tsuji), or they 
reflect a complicated time-scale competition (Allard & Homeier). 
The grains of different kind s have different distributions fs(a,z) 
in Marley, Ackerman & Lodders-model (i.e. Fe[s]-grains, H20[s] 
grains - different line styles in Fig.lH 3"* panel). Grains of different 
but homogeneous composition have the same size distribution in 
the Allard & Homeier and in the Tsuji-models at a particular height 
in the atmosphere. The dirty grains (i.e. a mixture of Fe[s]-SiO[s]- 
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Figure 3. Cloud material composition in volume fractions of total dust volume Vs/Vtot for prescribed model atmospheres of Teff =1800, 1400, 600K, log g = 
5.0 (Tsuji: Tc, = 1700K,Marley, Ackerman & Lodders: /sed = 2) resulting from different dust cloud models. 

Note: For Allard & Homeier VAUOs/^tot = (Vai203-cI + ^Al203-c2 + ^Al203-c3)/^tot and Vpe/Vtot = (Vpc-c + VFe-ci)/V,ot with c, cl, c2, c3 being different 
crystal structures treated in their equilibrium ansatz. For Tsuji no Teff = 600K model is available. 



Table 4. Maximum dust-to-gas ratios log {P(iust/Pgas)max> the maximum 
mean grain sizes log {a) [jim\ and its value in the upper cloud layers @ 
10'^ dyn/cm^ in different dust models for given (7, Psas, I'conv) profiles. 

Mai-ley, Allard& Helling& 

Teff Tsuji Ackerman Homeier Woitke 
& Lodders 

dust content 1800K -3.6 -3.2 -3.2 -3.4 

logPdust/Pgas 1400K -3.6 -3.2 -2.8 -3.4 

600K -5.8 (siiicmes) -4.5 -4.4 

-3.5 (H20) 



grain size 


1800K 


-2 


-1.4 .. 


- 1.6 





-2 


log {a) @ 


1400K 


-2 


-1.4 .. 


- 1.6 





-2 


10^ dyn/cm^ 


600K 




-1.6 .. 


-2 








maximum 


1800K 


-2 


+0.5 




+0.7 


+2.0 


grain size 


1400K 


-2 


1 




+0.5 


+2.5 


log <a>ma,K 


600K 


-2 


+0.5 




1 


-0.5 



Mg2Si04[s] etc.) in the Helling & Woitke-model are characterised 
by one mean grain size distribution f(a, z) at a particular height 
z in the atmosphere. The transition from (a>min to (a)niax across 
the cloud layer appears smoothly in the Helling & Woitke mod- 
els. The Marley, Ackerman & Lodders- models reach their <a>max 
abruptly at about the same latitude. Also the Allard & Homeier 
models show a sudden rise in mean grain size but at a different at- 
mospheric height. Another direct consequence of the different dust 
cloud modelling is that the values of (a) max can differ by a factor of 
100 (see also TableHll. 

3.2.3 Dust material composition 

The chemical composition of the cloud particles shows the largest 
variation between the different models. Figure[3]shows the material 



composition of the test models Teff = 1800, 1400, 600 K in volume 
fractions of the total dust volume, Vs/Vfot (definition Table [TJ, 
for the different cloud models. We only consider dust species which 
are important for the gas and dust opacity in the radiative transfer 
calculations in Sect. l3.3l (/' Table [Sjl. The calculation of the dust 
composition varies widely in the different models and splits into 
two classes: dust particles of homogeneous composition assuming 
equilibrium condensation (Tsuji; Marley, Ackerman & Lodders; 
Allard & Homeier) and dirty dust particles of heterogeneous com- 
position according to the kinetic treatment of growth and evapora- 
tion (Helling & Woitke). The chemical heterogeneity of the whole 
dust complex in the equilibrium models is reached by considering 
ensembles of pure Fe[s]-grains, Al203[s] -grains etc. 

The Tsuji and the Marley, Ackerman & Ackerman-models 
allow for the 3 and 5 major condensates, respectively, as dust 
opacity sources. The dominating low-temperature condensate is 
Mg2Si03[s] in the Tsuji-model and Mg2Si04[s] in the Marley, Ack- 
erman & Lodders-model. The Fe[s] fraction increases with de- 
creasing Teff in both models at intermediate temperatures. Also 
in the Helling & Woitke-model Mg2Si03[s] and Mg2Si04[s] are 
the dominant low temperature condensates in addition to Si02[s] 
and FeS[s]. Note that Si02[s] is never predicted by the equilibrium 
models. The Fe[s] content increases inward until it reaches, like 
in the other models, a pronounced maximum near the inner cloud 
edge. The dust at the cloud base is made of Al203[s] with pos- 
sible impurities of CaTiOsLs] (Teff=1800, 1400K) in the Helling 
& Woitke-model. The most refractory cloud condensate layer in 
Marley, Ackerman & Lodders-models is composed of corundum 
or Ca-Aluminates which serve as element sinks, hence, not de- 
picted in Fig. [3] The Allard & Homeier-models only partially agree 
with these results. The Marley, Ackerman & Lodders-model is the 
only model which includes H20[s] as possible condensate which 
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allows for a second, detached cloud layer above the already dis- 
cussed oxide-silicate cloud layer (compare Fig.|2j. 

Given the great diversity in grain composition with different 
model assumptions, we must conclude that the chemical composi- 
tion of the cloud particles in substellar atmospheres is still uncer- 
tain. 

3.2.4 Gas-phase chemistry results 

Figure |4] shows the number densities n [cm^^] for a selected num- 
ber of gas phase Si-, Mg-, A1-, Ca-, and Ti-bearing molecules. We 
additionally plot the most important H-molecules and the most im- 
portant C-bearing molecules (7"" and 8* panel). All models assume 
the gas-phase to be in chemical equilibrium {/^ Sect. |2. 11 1. All mod- 
els used the same well-mixed element abundances 6° at the inner 
boundary of the cloud mode {/^ Sect. l3.lt . Hence, different gas- 
phase number densities produced by the models are a consequence 
of the different treatment of dust formation which leads to different 
remaining element abundances ei in the gas phase. The compari- 
son of the remaining gas phase (after cloud formation) is needed 
to understand possible spectral trends in the later test cases of the 
complete (sub-)stellar atmosphere model (Sect. [33] l. 

3.2.4.1 H2 , Ht O, NH3 : The first test for differences in the chem- 
ical equilibrium gas-phase composition considers H2. Because of 
the continuous community interest we include H2O, and NH3 for 
its increasing spectral importance with decreasing Tefr in the sub- 
stellar regime. The H2 and H2O abundances are almost identical 
for all models. Figure |4] shows agreement also for the NH3 abun- 
dances except for the Allard & Homeier-model which predicts an 
overabundance of NH3 compared to the other models. 

3.2.4.2 CO, CH4: These molecules are only little effected by 
dust formation, since carbon solids are not considered in the mod- 
els under investigation. Hence, they are a good test for the general 
agreement of the gas-phase composition with respect to element 
abundances and material constants. However, the consumption of 
oxygen by the silicates and oxides does also affect the amount of 
gas-phase CO, hence indirectly also CH4, due to oxygen depletion. 
We observe that all models predict CH4 to be the major C-bearing 
molecule above a certain height in the atmosphere below which 
CO takes over. Beside this general agreement amongst the models, 
the CO number densities differ above the cloud layer, most likely 
resulting from different equilibrium constants for CO and CH4. 

3.2.4.3 TiO, Ti02: Ti02 is more abundant than TiO in all mod- 
els though the relative difference varies amongst the models. The 
models do not agree on the values of the TiO and Ti02 abundances. 
The Helling & Woitke model suggests the highest abundances for 
both molecules, the Tsuji-model suggests the lowest abundances. 

3.2.4.4 SiO, Si02: SiO is more abundant than Si02 in all mod- 
els though the relative difference varies widely amongst the mod- 
els. All models agree well for p > lO'dyn cm"^ which coincides 
with the pressure-level of the maximum dust content in this model 
(compare Fig.[2l left). The Tsuji-model again suggest the lowest 
molecular abundances, and the Helling & Woitke model suggest 
the highest abundances at lower pressures. 

3.2.4.5 MgH, MgOH: MgOH is more abundant than MgH in 
all models and the number densities agree well in the inner atmo- 
sphere for p > 10''dyn cm"^. The molecular abundances of MgH 



and MgOH fall into two groups with respect to the upper atmo- 
sphere: The Tsuji- and Helling & Woitke-models suggest a high 
number density. The Allard & Homeier- and Marley, Ackerman & 
Lodders-model predict the lowest abundances. 

3.2.4.6 Fe, FeH; AlH, AlOH: All models suggest that atomic 
Fe is more abundant than molecular FeH, and AlOH than AlH 
throughout the atmosphere and agree very well in the inner at- 
mosphere p > lO'dyn cm"^. All phase-equilibrium models yield 
very good agreements for all four molecules, and the Helling & 
Woitke suggests the highest Fe, FeH, AlH, and AlOH gas-phase 
abundances in the outer atmosphere layers. 

3.2.4.7 CaH, CaO, Ca(OH)2 : All models suggest that Ca(OH)2 
is more abundant than CaH for p < lO^dyncm^^. CaO has gen- 
erally a much lower abundance than these two molecules. The 
Helling & Woitke-model produces exceptionally high abundances 
of the Ca-molecules as result of the limited number of Ca-bearing 
solids i/' Sect. |2j2!4j. The Tsuji-model suggest the next highest 
abundances for these molecules. 

3.2.4.8 General: The general trend is that the phase-equilibrium 
models (Tsuji, Allard & Homeier, Marley, Ackerman & Lodders) 
produce lower gas-phase abundances of molecules containing dust- 
forming elements then the kinetic model (Helling & Woitke) in the 
upper atmosphere. However, the molecules that are not affected 
by the chemistry of dust formation (like CO, CH4, H2O) have 
very similar abundances in the different models. However, differ- 
ences for these molecules are indicative of the different oxygen- 
consumption caused by the differences in the dust cloud models, 
and of possible differences in the material quantities (Z' Sect. l2.lt . 
Those molecules containing rare element (like Al, Ti, Ca) are pre- 
dicted with very similar abundances in all phase-equilibrium mod- 
els. Remaining deviations for these molecules are likely due to a 
missing solid as element sink (Z' Ca Fig. |4j. The strongest devia- 
tions amongst these models occurs for molecules containing very 
abundant elements (Si, Mg). Since a large fraction of the Mg- and 
Si-bearing molecules contributes to the dust formation, the differ- 
ences in the dust models are imprinted in the remaining gas-phase 
abundances the strongest. All models agree on the gas-phase com- 
position below the cloud base (vertical black line. Fig. [2] left). 

3.3 Test case 2: global quantities given 

Two sets of stellar parameters were prescribed, 

T-dwarf Teff=1000K, log g = 5.0 
L - dwarf Teff = 1 800K, log g = 5 .0, 

for which complete model atmospheres were calculated including 
the solution of the radiative transfer. All models assume hydro- 
static equilibrium, gas-phase chemical equilibrium, and use mix- 
ing length theory for treating the convective energy transport. The 
dust cloud models are those described in Sect. [2] and all cloud ap- 
proaches assume spherical symmetric cloud particles. Table|5]sum- 
marised further details on the atmosphere codes. The cloud mod- 
ules are the model atmosphere component which is most different 
amongst the codes under consideration in this paper. 

3.4 Results test case 2 

We compare results for complete stellar atmosphere simulations re- 
garding the 
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Figure 4. Gas phase composition in a cloud forming atmosphere with Teff=1400K, logg=5.0 (Tsuji: Tcr = 1700K, Marley, Aclcerman & Lodders: /sej = 2) 
as resulting from different cloud model approaches. The vertical thin black line indicates the pressure at the cloud base where Pd/Pgas ~> {y Fig.|2] left). 



• atmosphere structure and cloud profile 

• spectral energy distribution 

• photometric fluxes and colours 



3.4.1 Atmosphere structure and cloud profiles 

Figure |5] shows the {T, p, v^.onv)-profiles and the cloud structures of 
the complete sub-stellar atmosphere simulations. 

The (r, p)-profiles ( P' panel) differ considerably in all parts of 
the atmosphere. Note that the Tsuji-model without dust opacity and 
the Tsuji-model with the highest Tc = 1900K are almost identical. 
The Dehn & Hauschildt -I- Helling & Woitke-model is the hottest 
throughout the entire atmosphere, where the Tsujix^^igooK. and the 
Tsuji-model without dust opacity, as well as the Allard & Homeier- 
model are the coolest at pressures below 10^ dyn cm"^. All models 
show the backwarming effect at T a 2000K for Teff = 1800K ex- 
cept the Allard & Homeier-model, the TsujiT„=i9ooK model, and the 
Tsuji-model without dust opacity. This backwarming is clearly as- 



sociated with occurrence of the cloud layer. It becomes stronger 
if the maximum amount of dust is situated at higher altitudes as 
the comparison of the different Tc -Tsuji-models demonstrates. The 
difference to the Allard & Homeier-model is understood by notic- 
ing that their model produces much less dust in the respective pres- 
sure regime where {a) (3''' and 4* row in Fig.|5j also Fig.[2]l 

The convective velocity v^onv explicitly enters all cloud model 
but the Tsuji-model (2"'' panel). The Allard & Homeier and the 
Dehn & Hauschildt -I- Helling & Woitke-model have identical vtonv 
since they use identical modules for treating the convective un- 
stable region below the Schwarschild limit but differ in the treat- 
ment of overshooting into the (classically) convective stable atmo- 
sphere at higher altitudes. The Tsuji-model and the Marley, Ack- 
erman & Lodders-model are comparable in v^a^y for the L-dwarf 
case Te(f = 1800K. Both suggest for the T-dwarf models a second 
convective layer which coincides with maximum Pdust/Pgas in these 
models (Fig. [5] right). 

The dust cloud structures are comparable in the sense that 
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Table 5. Summary of the model atmosphere codes used in Sect. |3.3l (s - soHds only, si - liquids and solids) 



authors 



code 
name 



element 
abundances 



number of number of number of 

elements gas-phase spec. dust species 



Tsuji 



AUard & Homeier 



Marley, Ackerman 
& Lodders 

Dehn & Hauschildt 



Anderse & Grevesse (1989) 
Allende Prieto et al.(2002) 

Settl- Grevesse, Noels & Sauval (1992) 
Phoenix Asplund, Grevesse & Sauval (2005) 



Lodders (2003) 
Grevesse, Noels & Sauval (1992) 



Drift- 



34 



84 



83 



40 



+ Helling & Woitke Phoenix 



83 3 as opacity source 

10 as element sinks 

680 43 as opacity source 

169 as element sinks 

- 2200 5 as opacity source 

~ 1700 as element sinks 

338 7 as opacity source 

7 as element sinks 



they appear in approximately the same pressure range except in 
the Tsuji-model where the small particles are homogeneously dis- 
tributed in the cloud until a certain critical temperature level Tent 
is reached. All models produce only the cloud layer of silicates 
and oxide^ (y Sect. 13.2.11 ). The results differ largely in the de- 
tails which confirms our results from the first part of our test cases 
(Sect. [T2l ). Different vertical cloud extension are suggested by dif- 
ferent simulations: the Allard & Homeier- model produces a ver- 
tically less extended cloud layer than the Marley, Ackerman & 
Lodders-model and the Dehn & Hauschildt + Helling& Woitke 
- model, which will have consequence for the emergent spec- 
trum of such an atmosphere. It appears that the differences in the 
cloud properties are amplified if the entire atmosphere problem is 
taken into consideration. The dust-to-gas ratio Pdust/Pgas (S"* row in 
Fig.|5} differs by 2 orders of magnitudes with the dust-opacity free 
Tsuji-model providing the upper limit and the Dehn & Hauschildt 
-I- Helling & Woitke-model, the Allard & Homeier and the 
TsujiT„=i7ooK - models the lower limit. Comparing this with the 
(r, p)-profiles (P' row) suggests that the higher the local tempera- 
tures at a certain pressure level for a given Teff=1800K, the smaller 

.Ackerman&Lodders ^ 

(Pdust/Pgas)Dchn&Hauschiidt+Heiiing& Woitke- TMs trend re-appears for 
Teff=1000K for which no Dehn & Hauschildt -I- Helling & Woitke- 
model is available. The Allard & Homeier-models are amongst the 
coolest of the (T, p)-profiles for a given Teff but always suggest a 

smaller pdust/Pgas- 

No trend appears in the numerical data for the mean parti- 
cle sizes (a). The models suggest the biggest particles to appear at 
the cloud base except the Tsuji-model which assumes a constant 
size in the entire cloud. Note that the particle sizes in the Allard & 
Homeier-models show a very steep distribution towards the max- 
imum size at the cloud base. It is apparent from Fig. |5] that the 
different dust cloud treatments produce different grain size distri- 
butions in the atmosphere resulting in differences in grain sizes up 
to 2 orders of magnitude at the cloud top where detectable spectral 
features would be produced (compare Marley, Ackerman & Lod- 
ders and Dehn & Hauschildt + Helling & Woitke). The Allard & 
Homeier-models suggest no particles in these cloud layers, hence, 
spectral dust features in these models might be more easily masked 
by molecular bands. 



^ The term silicate is loosely used among astronomers. In accordance with 
mineralogy, silicates are those solids containing Si-O groups. All other 
solids, like Al203[s], Ti02[s], CaTi03[s] etc. are oxides. 



3.4.2 Spectral energy distribution 

The spectral energy distributions (SEDs) between 0.5 ... 18/im 
calculated by the different model atmosphere codes employing dif- 
ferent cloud models are depicted in Fig. H (T^ff =I800K) and 
Fig. H] (Teff =1000K). The Marley, Ackerman & Lodders model 
for T^ff = 1000 K employed /^ed = 2 for consistency with the 
Teff = 1800 K case. Modeling by this group suggests that the spec- 
tra of early T dwarfs are better fit with larger values of fs^. The 
thick clouds resulting from the choice of /,ed = 2 (Fig. 5, bottom 
right) are responsible for the shallow absorption bands and red col- 
ors compared to the other groups for this case. In the Tsuji-case we 
plot two models, T;., = 1700K and = 1900K which demon- 
strate a very thin (T„ = 1900K) and an extended (T„ = 1700K) 
cloud layer. No two codes produce identical SEDs. Generally, 
Allard & Homeier and the Tsujixc,=i9ooK-models appear brighter 
than the other models between ~ 0.8 . . . ~ 1.5yum in the optical 
and near-IR (see also left panels Fig.[8]l. The Dehn & Hauschildt 
+ Helling & Woitke, the Marley, Ackerman & Lodders and the 
TsujiT(.r=i7ooK-models are the brightest of all models T^ff =1800K 
in the IR for A > 5//m. This result is not surprising because the 
Dehn & Hauschildt -I- Helling & Woitke, the Marley, Ackerman & 
Lodders and the TsujiTci=i700K-niodels contain the highest amount 
of small dust particles in the upper cloud layers (Sect. I3.4.1t and 
should therefore produce a redder atmosphere compared to a model 
without dust at comparable atmosphere pressures. In principle, the 
same analysis applies for the Teff =IOOOK-case, representing the 
T-dwarf regime within this study (Fig.|7]l. The cloud has moved al- 
ready considerably below t = I that both Tsuji-models appear very 
similar. Additionally, the Allard & Homeier-model suggests more 
spectral flux in several wavelength intervals than the Tsujixcr=i9ooK- 
model with an extended cloud layer Figure|7]also demonstrates an 
appreciable difference of the Tsuji-models around 10/jm which are 
sensitive to the application of the lOLA opacity band model. 

It is apparent from Fig. [8] that the Tsujixcr^isooK, the Mar- 
ley, Ackerman & Lodders-models and the Dehn & Hauschildt + 
Helling & Woitke-model produce much shallower absorption fea- 
tures in particular in the optical and near-IR then the TsujiTcr=i7ooK 
and the Allard & Homeier-models. Again, the reason lies in the 
differences in cloud modelling. This result appears surprising for 
the Dehn & Hauschildt + Helling & Woitke model since more el- 
ements remain in the gas phase due to incomplete dust formation, 
hence, the absorption features should be deeper A comparison with 
Fig. m (top panel) resolves this: All models showing shallow ab- 
sorption features are amongst the hottest (T,p) structures, hence, 
have low densities at given atmospheric temperature and there- 
fore lower opacity at that atmospheric height. Note that the (T, p)- 
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Figure 5. Test cases for complete atmospheric models for log g = 5.0, solar element abundance with Teff = 1800K (left) and Teff = lOOOK (right). 
Note: Different colours stand for different stellar atmospheres codes. Four models are plotted for the Tsuji-case (brown): long-short-dashed: Td =1700K 
(extended cloud), short-dashed: T^r =1800K, long-dashed: T^r =190()K (thin cloud), dotted: no dust opacity considered. Different line styles in log(a) indicate 
different homogeneous dust species in the Marley, Ackerman & Lodders-models. 



profiles are much more similar for Teff = lOOOK and so is the depth 
of the absorption features. Consequently, a spectral analysis relying 
on the depth of near- and IR spectral features, e.g. as gravity indi- 
cator, would underestimate the gravity in the case of the Tsuji and 
the Allard & Homeier-models compared to the Marley, Ackerman 
& Lodders and Dehn & Hauschildt + Helling & Woitke models in 
the L - dwarf regime. 

It appears that the treatment of the gas-phase opacity can ac- 
count for some differences in the synthetic spectral energy distri- 
bution, here in particular the treatment by the JOLA band method 
in the Tsuji-models vs. the more frequency sensitive methods used 
in all other models. However, the completeness of the molecular 
line lists has only minor effects on our photometry results given the 
large influence of the dust modelling demonstrated here. 



3.4.3 Photometric fluxes 

We wish to compare our simulations also in terms of photomet- 
ric fluxes (Figs.[8]|9lll0|and Table [6j. We have chosen to demon- 
strate our comparison for four filter systems covering the near- 



IR and the IR: the JHK-2MASS photometric systenfl the WF- 
CAM UKIRT filterfl the IRAC Spitzer photometric band;]!, and 
the VISIR VLT systenfl Delfosse et al. (2000) show for their em- 
pirical mass-luminosity relation that model atmosphere results are 
more reliable in the near-lR (JHK) than at lower wavelength. How- 
ever, we include the Z and Y (+JHK) from WFCAM UKIRT filter 
system for comparison. Note that Carpenter (2001) provides trans- 
formation formula for the 2MASS colours into a number of differ- 
ent photometric systems (also Hewett et al. 2006). 

The Allard & Homeier-model exhibits the highest J, H and 
Y photometric fluxes in the Teff = 1800K-case, while both 
Tsuji-models bracket the compared models for the IRAC pho- 
tometric fluxes (except Band 4). The Allard & Homeier, and 
the Tsujix„=i9ooK-models have the largest Y fluxes in the Teff = 
lOOOK-case, and the Tsujix„=i7ooK-model has the largest flux in Z. 
The Marley, Ackerman & Lodders-model and the TsujiT„=i9ooK- 

* http://web.ipac.caltech.edu/staff/waw/2mass/opt_ cal/index.html 

^ http://www.ukidss.org/technical/technical.html 

^ http://ssc.spitzer.caltech.edu/irac/spectral_ response.html 

' http://www.eso.org/sci/facilities/paranal/instruments/visir/inst/index.html 
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Figure 6. Synthetic spectra for Teff =1800K, logg=5.0 and solai' metalicity. Two spectra are plotted for the Tsuji-model: Tc, 
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= 1700K (brown; extended cloud) 



model suggest the largest flux in Z for T^ff 
7UKIRT ^ logFj^^^^^, logFj^™'^ a 



= 1800K. Note 
jQg^2MASS^ and 



that logF} 

log F^''"*^ < log F^^^^ for all model approaches. The model re- 
sults diff^er the most in the UKIRT ZYJ bands in both test cases. 

Interestingly, the Tsuji models suggest the faintest fluxes in 
all IRAC band for T^g = lOOOK while the Maiiey, Ackerman & 
Lodders-models result in the largest fluxes in these wavelength 
bands. The maximum difference amongst the models in photomet- 
ric fluxes (An,ax[logFc], last columns Table [6} are larger in the T- 
dwarf test case than in the L-dwarf test case. The maximum dif- 
ferences occur in the ZYJH bands, and in the VISIR SIV, SiC, and 
Nell bands for the T-dwarf test case. 

These synthetic photometry allows to suggest an error margin 
due to spread in the model results for apparent magnitudes (defi- 
nition see Table [TJ. The photometric flux differences relate to un- 
certainties in apparent magnitudes between 0.25 < Am < 0.875 
for the L-dwarf test case. The uncertainty in apparent magnitudes 
increases considerably for the T-dwarf test case due to strong dif- 
ferences in the Y band: 0. 1 < Am < 1.375 . 



Figures |9] and [TOl demonstrate how the photometric fluxes of 
the model atmosphere codes translate into colours (definition see 
Table [TJ. For this, we have normalised the photometric fluxes to 
the corresponding photometric fluxes for Vega (log F^q in Tables [T] 
and|6j. The largest differences occur for the L-dwarf model (Teff = 
1800K) in J-[4.5], Ks-[4.5], K-[4.5], and the [3.5]-[4.5]. For the 
T-dwarf model (Teff = lOOOK) the maximal differences occur in 
Ks-[3.6], [3.5]-[4.5], and Ks-[4.5]. 

If we assume the synthetic colours are correct, we can in prin- 
ciple use them to infer the spectral type of an object from an ob- 
served (colour, SpT)-diagra n0. Ideally, the synthetic colours de- 
rived in Fig. [9] for Teff=1800K should result in a L-dwarf and the 
colours for Teff=1000K should suggest a T-dwarf. For this exercise, 
the synthetic colours (mf - m2) in Fig.[9]are compared with the ob- 
served IRAC colours in Patten et al. (2006), and the related spectral 
type is picked from their Fig. 10. However, in Patten et al. (2006) 
the [3.6]-[4.5]-SpT relation appears ambiguous for M7 . . . L8 since 



SpT = Spectral Type 
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Figure 7. Synthetic spectra forTeff =1000K, logg=5.0 and (initial) solar metalicity. Two spectra are plotted for the Tsuji-model: Tc, =1700K (brown; extended 
cloud) and T^r =1900K (orange; thin cloud). 



[3.6]-[4.5] remain approximately constant for these spectral types. 
Our synthetic error margin is slightly larger than the [3.6]-[4.5] 
scatter in Patten et al. (2006). We encounter a similar challenge in 
relating our synthetic mean value for [4.5]-[5.8] to a possible spec- 
tral type suggesting MO . . . L7 in the worst case, but the mean value 
for [4.5]-[5.8] resuhs in SpT= M8 . . . L6. Patten et al. (2006) ar- 
gue that the scatter in their observed [5.8]-[8.0]-SpT plot might be 
due to the H2O vs. CH4 absorption at 7.7/jm. Our synthetic mean 
value would then suggest with these observed data an interval of 
M9 . . . T5 inside the synthetic error bars, while the mean value for 
[5.8]-[8.0] narrows the SpT range to L0...T4. About the same 
conservative spectral type range is found for J-[4.5], but can be 
narrowed to L3 . . . L5 for the J-[4.5]-mean value. For our T-dwarf 
test case (Teff = lOOOK) the biggest uncertainties in comparison 
with the Patten et al (2006) data occur for [5.8]-[8.0], Ks-[3.6] and 
J-[4.5]. However, the SpT-colour relation in Patten et al. (2006) is 
much narrower for T-dwarf s, and hence, our synthetic mean colours 
do suggest much narrower SpT ranges than for our L-dwarf test 
case (compare Table|7j. 



Interestingly, the two extreme Tsuji-models (T^., = 1700K and 
T„ = 1900K) bracket the Y-J and the Z-J UKIRT colours in FigfTOl 
We compare our synthetic Y-J and J-H UKRIT colours to Hewett et 
al. (2006), and we reproduce the spectral class of our test case mod- 
els better than for the near-IR colours (see Table [TJ. A comparison 
with Lodieu et al.(2007a) demonstrates that synthetic Z-J colours 
falls well in their sequence of substellar objects with decreasing 
mass. However, a proper reproduction of a T-dwarf spectrum would 
demand an adjustment of cloud parameters like T„it, f^i and pos- 
sibly also the mixing efficiency as described in Sections |2.2. 1.21 
12.2.2.21 l2.2.3.4l and l2.2.4.4l compared to the parameter used in this 
comparison study which are more suitable for L-dwarf model at- 
mospheres (see e.g. Tsuji 2005, Gushing et al. 2008, Stephens et al. 
2008). 



3.4.4 General 

We have evaluated our model results in the most conservative 
way and despite the differences amongst the simulations in mod- 
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Table 6. Photometric fluxes log [erg cm^^s"' A"'] with c indicating the JHK-2MASS system, the WFCAM UKIRT system, the VISIR-bands, and the IRAC 
Spitzer bands for the models depicted in Fig. [8] The last column contains the maximum differences in log amongst the models: Aniax[log ^c] = [log fc]'"" ~ 
[logfc]"""- Note that Tsuji's thick-cloud case (T^r = 1700K; each 1st row) is used to calculate A„,ax[logfc] for the L-dwarf test case (Teff = 1800K), and 
Tsuji's thin-cloud case (Tj, = 190()K; each 2nd row) for the T-dwarf test case ( Teff = lOOOK). We also list the photometric fluxes for the HST Vega spectrum 
of Bohlin & Gilliland (2004; Wan'en 2008, priv. com.) which we use as zero points to calculate the synthetic colours (y TablefT] Figs. |9] [TO). 
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elling clouds, they suggest a general agreement in synthetic colours 
(Figs.|9l not . The mean values do provide a good guidance taking 
into account the large diversity regarding the characteristic cloud 
quantities like the mean particle sizes and material composition. Of 
course, uncertainties increase if the error margins given in Table |7] 



are applied. The values of the spectral types suggested in Table [7] 
suffers also the scatter and possible ambiguity contained in the ob- 
served data like e.g. in the IRAC colours [3.6]-[4.5] or [5.8]-[8.0] 
as discussed by Patten et al. (2006) which adds to the model inher- 
ent uncertainties. A better picture appears if colours are correlated. 
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Figure 8. Synthetic spectra for Teff =1800K (top) and Teff =1000K (bottom) with logg=5.0 and solar metalicity for a spectral resolution R=200. Two spectra 
are plotted for the Tsuji dust model: Td =1700K (brown; extended cloud) and Tc, =1900K (orange; thin cloud). Photometric fluxes (symbols) are plotted 
for the JHK-2MASS-system (left panels, ■) and the IRAC Spitzer-bands (right panels, A) at the band center frequency AA/2. The photometric fluxes ai'e 
summarised in Table[6] 



A comparison with colour-colour plots (e.g. Lodieu et al. 2007b 
for UKIRT) shows that the different simulations reproduce well the 
spectral classes of the test cases even in the most uncertain ZYJ 
colours. Nevertheless, we refrain from the exercise of back-tracing 
the Tcff-values from our synthetic colours since this would clearly 
potentiate uncertainties because the atmosphere simulations used 
in publications of interest did not treat the presence of dust at all 
(in Luhman 1999) or used very simplistic representations of dust 
as opacity source (in Golomowski et al. 2004). 



4 DISCUSSION 

4.1 The challenge of phase-transition modelling 

An essential part of modelling clouds in substellar atmospheres 
is the description of condensation as a phase-transition gas - 
solid/liquid. Two modelling approaches were used in the simula- 
tions compared in this paper: the kinetic approach and the phase- 
equilibrium approach. Condensation occurs in the kinetic approach 
when a gas species is supersaturated with respect to its equilib- 
rium concentrate at given pressure and temperature. Homogeneous 
nucleation of a supersaturated species yields macro-molecules, 
molecular clusters and eventually nanometer solids. The clusters 
and nanometer solids could become seeds for heterogeneous nu- 
cleation by surface reactions that requires lower activation energies 
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Table 7. Mean synthetic colours ((mi - m2)n,eim = Zf ("^l ~ t>'2)l/L; L - to- 
tal number of models) and error margins (((mi - m2)max -{mi- 'M2)min)/2 - 
maximum colour difference) derived from Figsl9l [Tol The HST Vega spec- 
trum of Bohlin & Gilliland (2004; Warren 2008, priv. com.) is used as zero 
point. Listed are also the spectral types suggested by the mean synthetic 
colours alone (each 1st row) and including the synthetic error margin (each 
2nd row). 
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Figure 9. Synthetic colour-colour diagrams for 2MASS, for Teff =1800K 
(□ and cross) and Tj-ff =1000K (A and star) with logg=5.0 and solar metal- 
icity. The symbol-colour code is the same like in Fig.|6] Tsuji T„ =1700K 
(brown, thick cloud) and Tc, =1900K (orange, thin cloud), Marley, Ack- 
erman & Lodders (blue), Allard & Homeier (green), Dehn& Hauschildt + 
Helling & Woitke (magenta). The star symbols indicate the colours evalu- 
ated for the WFCAM UKIRT filter system. 



than homogeneous nucleation of the same species directly out of 
the gas phase. However, the first condensate will form by homoge- 
neous nucleation. A second condensate can form at a lower tem- 
perature by heterogeneous nucleation on seeds made of the first 
condensate or by homogeneous nucleation following supersatura- 
tion of the gas. In the phase-equilibrium approach, thermodynamic 
equilibrium is adopted where the Gibbs energy difference between 
reagents and products equals zero at fixed values of temperature, 
pressure and element composition. A series of these calculable 
states when ordered as a function of decreasing temperature can be 
viewed as a time sequence of fractional condensation that predicts 
the stepped appearance of crystallographically ordered, chemically 
stoichiometric solids, i.e. minerals. Equilibrium condensation mod- 
els thus predict a sequence of minerals of systematically different 
compositions. Yet, the presence of a particular mineral or mineral 
assemblage is no proof of equilibrium condensation. 

Laboratory condensation experiments on silicate vapours 
found that dissipative structures (i.e. metastable states) ap- 
pear as highly disordered, amorphous solids with unique non- 
stoichiometric compositions. That is, they have unique Metal-oxide 



to Si02 ratios (IVI: IVIg, Fe, Ca, or Al, and combinations thereof) that 
match deep metastable eutectic compositions in equilibrium phase 
diagrams (Nuth et al., 1998, 2000; Rietmeijer et al., 1999, 2008). 
These compositions are intermediate between equilibrium min- 
eral compositions, e.g. serpentine dehydroxylate, ]VIg3Si207. With 
time it would breaks down to forsterite and enstatite (]VIg2Si04 
+ IVlgSiOs) (Rietmeijer et al., 2002) that both would form at dif- 
ferent temperatures during equilibrium condensation. Equilibrium 
condensation is predictable but so again is extreme non-equilibrium 
condensation of deep metastable eutectic condensates that are more 
reactive than equilibrium minerals. With time, post-condensation 
thermal annealing (i.e. ageing) of non-equilibrium condensates will 
also lead to thermodynamic equilibrium minerals. Ageing is de- 
termined by the prevailing time-temperature regime (Hallenbeck 
et al., 2000) and condensate morphology, i.e. aggregates or dust 
clumps (Rietmeijer et al., 1986, 2002). The ageing process will 
yield minerals at temperatures below their equilibrium condensa- 
tion temperature. Such processes would require the cloud particles 
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Figure 10. Synthetic colour-colour diagrams only for WFCAM UKIRT fil- 
ters for Teff =1800K (cross) and Teff =1000K (star) with logg=5.0 and 
solar metalicity. The symbol-colour code is the same like in Fig. [6] Tsuji 
Tcr =1700K (brown, thick cloud) and T^, =1900K (orange, thin cloud), 
Marley, Ackerman & Lodders (blue), Allard & Homeier (green), Dehn& 
Hauschildt -t- Helling & Woitke (magenta). 

to remain in a certain thermodynamical state long enough, a situa- 
tion possibly occurring to dust trapped in-between convection cells. 

Toppani et al. (2006) demonstrate that their (Mg,Ca,Al,Si)- 
oxide vapour condenses to complex hydrates carbonates in a CO2- 
H20-rich gas, and they conclude that this condensation proceeds 
near-equilibrium. Their condensation experiments at moderate gas 
temperate and low total pressures (1000-1285K, 0.004 bar) yield 
many of the expected equilibrium condensates in crystalline form. 
However, both Toppani et al. (2006) and Rietmeijer et al. (2008) 
conclude that the mineralogy of such condensed material can not 
be understood without taking into account the influence of kinetics. 



5 CONCLUSION 

Clouds in the atmospheres of brown dwarfs and gas-giant planets 
determine their spectral appearance and influence their evolution by 
altering the atmospheric thermal structure. The challenge of mod- 
elling cloud formation has been approached from very different per- 
spectives over the past years which leads to the question: Do these 
models yield the same results and how much do they differ in pre- 
dicted observational quantities? Five models are compared in this 
paper to address these questions. All models emphasise the chem- 
istry of cloud formation. Considering clouds as the result of a phase 
transition process (gas to solid/liquid), the models assuming phase- 
equilibrium describe the end-state of the phase-transition process, 
whereas kinetic models describe the initial state of the cloud for- 
mation process from a chemical point of view. Which viewpoint is 
most correct ultimately depends upon the timescales for the various 
relevant atmospheric processes. 

The dust cloud models predict generally comparable cloud 
structures despite the different approaches, although the results dif- 
fer substantially in detail. Opacity relevant quantities like grain 
size, amount of dust, dust- and gas-phase composition vary be- 
tween the various approaches. Most cloud models agree that small 
grains composed mainly of silicates (MgSi02[s]/Mg2Si04[s]) pop- 
ulate the upper cloud layers, whereas iron (Fe[s]) is a major com- 



ponent of the large grains at the cloud base. The cloud models 
agree on the gas-phase composition in the inner atmosphere only 
which is too warm for condensed phases. All models predict phase- 
equilibrium here, though the different models describe the evapora- 
tion at different levels of detail. Above the cloud, more molecules 
remain in the gas-phase if cloud formation is treated in phase-non- 
equilibrium compared to results from phase-equilibrium models. 
The dilferent results that arise from differences in cloud modelling 
are amplified if the entire atmosphere problem is solved, including 
radiative and convective energy transport. The reason is the strong 
feedback of the clouds on the (r, /7)-structure due to the clouds' 
strong opacity and its high efficiency in depleting the gas of the 
atmosphere. 

Viewing their spectral appearance, the results of the cloud 
model atmosphere codes appear to fall into two categories: 

- The high-altitude chmd models or extended (Tsujixcr=i7ooKj Mar- 
ley, Ackerman & Lodders, Dehn & Hauschildt -1- Helling & 
Woitke) where the dust-to-gas ratio peaks at high altitudes though 
at different absolute levels. In these models, small grains ({a) = 
10"* . . . lO^'^fim) are still present well above the maximum of 
Pd/Paas. hence the gas-phase absorption is less deep for A > lyum. 

- The low-altitude cloud models or thin (Tsujixc-igooK. Allard & 
Homeier) where the dust-to-gas ratio maximum sits further inside 
the atmosphere and no grains populate higher atmospheric layers 
above the maximum of Pd/Pgas- Consequently, the gas-phase ab- 
sorption features are much deeper in the mid-IR and IR part of the 
spectrum. Consequently, the low-altitude cloud / thin model atmo- 
spheres appear bluer then the high-altitude cloud / extended model 
atmospheres. 

Comparing synthetic photometric fluxes and colours from dif- 
ferent model atmosphere codes illustrates the current range of un- 
certainty, or error bar, for theoretical predictions. These error bars 
are worst cases. They are derived from a group of different mod- 
els and not from only one particular family of models. In the 
most conservative case, the maximum differences in photomet- 
ric fluxes, A,^ax[log amongst the models are between 1% and 
30%. An,im[logFi] increases to 50% for the T-dwarf test case in 
the WFCAM UKIRT wavelength intervals. This translates into an 
uncertainty in apparent magnitudes for the L-dwarf test case of 
0.25 < Am < 0.875. Ani increases to 1.375 for the T-dwarf test 
case in the WFCAM UKIRT filter system. 

We conclude that every comparison with observations should 
ideally involve models from different groups. This would allow the 
determination of a synthetic error bar in determinations of funda- 
mental quantities like Teff, log g, and metallicity. Ultimately com- 
parison of models to objects with independently constrained prop- 
erties (from orbital motion and bolometric luminosity, for example) 
will elucidate the modelling approaches that most accurately cap- 
ture the relevant physics. Other possibilities for tests include ob- 
jects in young stellar clusters which have well constrained ages and 
metallicities. 

Future works on cloud formation need to seek more support 
in laboratory astrophysics. Hydrodynamic modelling of ultra-cool 
atmospheres will provide the opportunity to study the dynamic pro- 
cesses of cloud formation. They need to including a consistent de- 
scription of the chemical formation processes and simultaneously 
address the challenge of a turbulent fluid field. 
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